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Abstract 

The bosonic su(n) Hubbard model was recently introduced. The model was 
shown to be integrable in one dimension by exhibiting the infinite set of conserved 
quantities. I derive the i?-matrix and use it to show that the conserved charges 
commute among themselves. This new matrix is a non-additive solution of the 
Yang-Baxter equation. Some properties of this matrix are derived. 



PACS numbers: 75.10.-b, 75.10. Jm, 75.10.Lp 

Key words: Hubbard model, su(n) spin-chain, integrability 



October 6 th 1997 
LAVAL-PHY-25/97 



*Work supported by NSERC (Canada) and FCAR (Quebec), 
^email address: zmaassar@phy.ulaval.ca 



1 Introduction 



The two-dimensional Hubbard model was introduced to describe the effects of correlation 
for (^-electrons in transition metals 0. It was then shown to be relevant to the study of 
high-T c superconductivity of cuprate compounds. 

In one dimension the model is integrable || |J. The integrability framework of 
the model is the quantum inverse scattering method [H. However, despite sharing many 
properties with discrete quantum integrable models, the model has a peculiar integrable 
structure which defines a class of its own. 

In seeking to generalize the Hubbard model in any dimension, it was therefore natural 
to look for a one-dimensional generalization which is integrable. An n-state generalized 
model which contains the usual su(2) model was recently introduced in ||. This su(n) 
Hubbard model was shown to possess an infinite set of conserved charges and to have 
an extended su(n) symmetry. The model is built by coupling two copies of the recently 
discovered su(n) XX 'free-fermions' model 0. For n = 2 a fermionic formulation exists, 
but for n > 2 finding an analogous framework is a tantalizing problem. 

In this work I derive the i?-matrix of the model; this provides a direct proof of the 
commutation of the conserved charges among themselves. Section two gives the definition 
of the bosonic Hamiltonian and the transfer matrix. The i?-matrix intertwining the 
monodromy matrices is derived in section three. In section four some properties of this 
new matrix are given. I conclude with some remarks and outline some outstanding issues. 



2 The model 

Let E a @ be the nx n matrix with a one at row a and column j3 and zeros otherwise. The 
su(n) Hubbard Hamiltonian on a ring then reads ||: 

h 2 = E^+i+EC + ^E^ (!) 

i i i 

= E E (^rE™ + x^E™E™ + (E^ E')) + U^( Pi + + ^) 

i a<n i 

where p = J2 a <nE aa — (n — l)E nn , and primed and unprimed quantities correspond 
to two commuting copies of the E matrices. The Hamiltonians h and h are su(n) XX 
Hamiltonians [Tjj. The complex free parameter x is a deformation inherited from the XX 
model. The Hamiltonian H2 is defined in one dimension but can be evidently defined on 
any lattice; integrability is lost however. 

For n = 2 and x = 1, and using Pauli matrices, the Hamiltonian is just the integrable 
bosonic version of the usual Hubbard Hamiltonian || : 

h? ] = |E(«i + °Mn) + (* -> <0 + UY,°t°i (2) 

/ i i 

The Hamiltonians can be written simply in terms of su(n) hermitian traceless matrices. 
For |x| = 1 the Hamiltonians are hermitian. 

The transfer matrix is the generator of the infinite set of conserved quantities. Its 
construction was given in 0. We recall it here. Consider first the i?-matrix of the su(n) 
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XX model [fl: 



R(\) = a(X) [E nn (g)E nn + E<ia ® E ° P ] 

a,0<n 

+6(A) ^ (xE nn ® £ Qa + aT 1 ^ ® £ r ' 



+c(A) (E na ® i? on 4- i? Qn ® E na ) (3) 



where a(A) = cos(A), 6 = sin(A) and c(A) = 1. The functions a, b and c satisfy the 
'free-fermion' condition: a 2 + b 2 = c 2 . For this set of parameters, a Jordan- Wigner 
transformation turns the U — Hamiltonian density for su(2) into a fermionic expression 
for free fermions hopping on the lattice. 
Consider also the matrix 

I Q (h) = cosh(^) Id + sinh(^) Co C' = exp ^ C C'A (4) 

where C = J2 a <n E aa ~ E nn . We stress that C turns out to be the fundamental matrix, 
not the su(n) generator p. We have p + ^j^Id = |C, for n > 2. The parameter h is 
related to the spectral parameter A by 

n 2 U 

sinh(2/i) = — sin(2A) (5) 

One chooses for h(X) the principal branch which vanishes for vanishing A or U. Then for 
U = the monodromy matrix becomes a tensor product of two uncoupled XX models. 
The Lax operator at site i is given by: 

L 0l (X) = I (h) R 0l (X) R' 0l (X) I (h) (6) 

and the monodromy matrix is a product of Lax operators, T(X) = Lom(A)...Loi(A), where 
M is the number of sites on the chain. The transfer matrix is the trace of the monodromy 
matrix over the auxiliary space 0: r(A) = Tr [(L 0M ...L 01 ) (A)]. One possible set of 
conserved quantities is given by 

Hp+1 = (fhiR) (7) 



dXP 



A=0 



The proof that H2 commutes with r(A) was given in ||. The derivative of the matrix / 
gives the coupling term appearing in (0). Note that the definition involving a logarithm 
has two benefits. Besides giving the most local operators, it further disentangles the two 
copies. 



3 The i?-matrix 

We derive the i?-matrix intertwining two monodromy matrices at different spectral pa- 
rameters. To this end we generalize the algebraic method of the Decorated Star Triangle 
Equation introduced by Shastry [TjJ . 
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The XX i?-matrix satisfies the regularity property R(0) = Id, the unitarity condition 
R(X)R(—X) = Id cos 2 A and the Yang-Baxter equation 

R 12 (X - fi)R 13 (X)R 23 (/j) = RMR 23 (\)R 12 (\ - /j) (8) 

where R = PR and P is the permutation operator on the tensor product of two n- 
dimensional spaces. It is easy to verify that it also satisfies a decorated Yang-Baxter 
equation 

Ri 2 (X + fi)C 1 R 13 (X)RM = RMR 23 (\)C 2 R 12 (\ + f i) (9) 
We now look for the i?-matrix intertwining two L-matrices: 

R(X 1 , A 2 ) L(Ai) ® L(A 2 ) = L(A 2 ) <g> L(A X ) fl(A x , A 2 ) (10) 

The su(2) case lead us to consider the following Ansatz [jT0| : 

^(Ai,A 2 ) = /i 2 (/i 2 )/ 3 4(/ii)(«^i3(Ai-A 2 ) J R 24 (A 1 -A 2 ) 

+/3 J R 13 (A 1 + A 2 )C 1 A 24 (A 1 + A 2 )C 2 ) h 2 {-h 1 )h i {-h 2 ) (11) 

The .R-matrix acts on the product of four auxiliary spaces labeled from 1 to 4, and a, (3 
are to be determined. One then requires relation ([TOj) to be satisfied and uses (||) and (|9]) 
to derive the following equation: 

(a J R 13 (A 1 - A 2 )i? 24 (A 1 - A 2 ) + P^Rn^ + A 2 )C 4 i? 24 (A 1 + A 2 )) Ji 2 (2/n)/ 34 (2/i 2 ) = 

/i 2 (2/i 2 )I 3 4(2/ii) (a i?i 3 (Ai - A 2 ) J R 24 (A 1 - A 2 ) +/3R 13 {X 1 + A 2 )C lJ R 24 (Ai + A 2 )C 2 ) 
Expanding the exponentials and taking into account all the terms yield only two equations: 

— = — tanh(/i! + h 2 ) , — = — tanh(/i! — h 2 ) (12) 
a B a A 

where a = cos(Ai — A 2 ), b = sin(Ai — A 2 ), A = cos(Ai + A 2 ) and B = sin(Ai + A 2 ). The 
compatibility equation 

tan(Ai - A 2 ) _ tanh(/i x - h 2 ) 

tan(Ai + A 2 ) tanh(/ii + h 2 ) 
is satisfied provided equation (0) is satisfied for the pairs (Ai,/ti) and (X 2 ,h 2 ). One can 
then pull out a = a(X\, A 2 ) which appears as an arbitrary normalization of the i?-matrix, 
to obtain: 

R{X U A 2 ) = a(Ai, X 2 )I 12 (h 2 )I u (h) { Rn{X x - A 2 )i? 24 (A 1 - A 2 ) + sm [ Xl - ^ 

\ sm(Ai + A 2 ) 

x tanh(/n + h 2 )Ri 3 (X 1 + A 2 )C lJ R 24 (Ai + X 2 )C 2 ) J 12 (-/ii)J 34 (-/i 2 ) (14) 

The monodromy matrix being a tensor product of M copies of L matrices, one has 

R(Xx,X 2 ) T(Ai) ® r(A 2 ) = T(A 2 ) ® T(Ai) £(A ls A 2 ) (15) 

Taking the trace over the auxiliary spaces and using the cyclicity property of the trace 
one obtains [r(Ax), r(A 2 )] = 0. We have thus proven that all the conserved charges H p 
mutually commute. 

Note that this proof is rigorous and valid for all values ofn, and for arbitrary values of 
the complex parameter x. It only involves the algebraic properties of the operators appear- 
ing in the various matrices, not the specific n-dependent representation. The equations 
(H) and (P) are the only equations of this type needed for the proof. 
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4 Properties of the R matrix 



I now give some properties of the i?-matrix. At XJ = the two XX models decouple 
and h(X, U) = 0. Expression flT4| ) indeed decouples as a tensor product of two su(n) XX 
^-matrices. 

The matrix also satisfies the regularity property 

R(Xi, Ai) = a(\x, Xx) Id (16) 

and the unitarity property: 

R{X 1 ,X 2 )R{X 2 ,X 1 ) = a 2 (A 1 ,A 2 )cos 2 (A 1 -A 2 ) 

x (cos 2 (Ai - A 2 ) - cos 2 (Ai + A 2 ) tanh 2 (/i 1 - h 2 )) Id (17) 

The derivation of the last property is straightforward and involves algebraic relation be- 
tween the building blocks of the su(n) XX R- matrix. 

One can invoke the associativity of the algebra of L matrices, which ultimately re- 
duces to the associativity of usual matrix multiplication, to conclude that the intertwiner 
satisfies a Yang-Baxter relation of its own. The two ways of permuting a product of three 
L-matrices imply 

ft L(Ai) ® L(A 2 ) ® L(A 3 ) IT 1 = L(A X ) <g> L(A 2 ) <g> L(A 3 ) 
ft = (.Ri 2 (A 2 , A 3 )^ 23 (Ai, A 3 )-Ri 2 (Ai, A 2 )^) -R 23 (Ai, A 2 ).Ri 2 (Ai, A 3 )-R 23 (A 2 , A 3 ) (18) 

I am unaware of the existence of an equivalent of the Schur lemma for the algebra of L- 
matrices. This would allow to conclude that II oc Id. Once proportionality is established, 
the regularity property ensures that the proportionality constant is one. We can argue 
that II = Id holds because it has been explicitly verified for n = 2 ||, and because the 
building blocks of the matrix satisfy algebraic relations which are independent of n. Thus 
the i?-matrix satisfies the Yang-Baxter equation: 

i?i 2 (A 2 , A 3 ).R 23 (Ai, A 3 )i?i 2 (Ai, A 2 ) = .R 23 (Ai, A 2 )i?i 2 (Ai, A 3 )i? 23 (A 2 , A 3 ) (19) 

where A and h are related through @. 

Using the explicit expression (|14|), it should be possible and it is instructive to try to 
check that the above equation is satisfied for any value of n. The factors / drop out and 
half of the terms on both sides of the YBE compensate each other because relations (||]) 
and (^) hold. The eight remaining terms involve highly non-trivial relations. Each term 
is a product of six /^-matrices, three for every copy, in the ordering dictated by the YBE. 
The C factors can be dropped by changing the arguments of the i?-matrices appropriately. 
However the arguments do not allow the use of (|8], ||) because the middle argument is not 
the sum of the extreme ones. One is forced to expand all the products on a basis and to 
regroup terms and check that the resulting trigonometric constraints are satisfied. I have 
not verified whether all these relations hold. Although specific particularities pertaining 
to the XX matrix are needed, I stress again that the proof is algebraic. In this respect 
the proof of || for n = 2, although following a different approach, should generalize in a 
straightforward way to any value of n. 
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5 Conclusion 



We have shown that all the conserved charges of the su(n) Hubbard model mutually 
commute by exhibiting the intertwining matrix. This matrix is the su(n) generalization 
of the su(2) one obtained in ||. Some properties where then derived. A notable feature 
of the matrix is its non-additivity property; the A dependence cannot be reduced to a 
difference (Ai — A 2 ). 

One can now start diagonalizing the Hamiltonians by the method of the algebraic 
Bethe Ansatz. Preliminary results suggest an interesting structure for the Bethe eigen- 
states p. 

One can also consider the extension of the su(n) model to other algebras. The algebraic 
underlying structure of the su{n) XX model should admit generalizations [|7|, [|. 

Acknowledgement: I thank P. Mathieu for bringing to my attention reference flTC| ] 
and for continued support. 

Note: While this work was being written, Martins exhibited a gauge-transformed 
version of the foregoing R(\i, A2) matrix. The derivation is also based on a generalization 



of Shastry's method. However, the proof in [0 was carried out for n = 3, and 'extensive 



checks' were made for n = 3,4. The expression of R for all values of n is left as a 
(correct) conjecture. This reference also used unnecessarily complicated versions of the 
Yang-Baxter equations. 
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